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Abstract 

' We perform the periodic homogenization (i. e. e 0) of the non-stationary Stokes-Nernst-Planck- 

(-H \ Poisson system using two-scale convergence, where e is a suitable scale parameter. The objective 

, is to investigate the influence of difFerent boundary conditions and variable choices of scalings 

in e of the microscopic system of partial differential equations on the structure of the (upscaled) 
limit model equations. Due to the specific nonlinear coupling of the underlying equations, special 
^ \ attention has to be paid when passing to the limit in the electrostatic drift term. As a direct result 

, of the homogenization procedure, various classes of upscaled model equations are obtained. 
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1. Introduction 

5t 1 This paper deals with with the periodic homogenization of a non-stationary Stokes-Nernst- 

Planck-Poisson- type system (SNPP). The real- world applications that fit to this context include 
areas of colloid chemistry, electro-hydrodynamics and semiconductor devices. Our interest lies 
in the theoretical understanding of colloid enhanced contaminant transport in the soil. Colloidal 
particles are under consideration for quite a long time since they are very important in multiple 
applications ranging from waste water treatment, food industry, to printing, etc. The monograph of 
van de Ven ^3(1] and the books by Elimelech [ll | and Hunter [isl yield a well founded description of 
colloidal particles and their properties. However, the difFerent processes determining the dynamics 
of colloids within a heterogenous porous medium are not yet completely understood. Therefore, 
the mathematically founded forecast of contaminant transport within soils is still very difficult, as 



it is strongly influenced by the movement and distribution of colloidal particles (cf. e.g. [29|). 
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Using mathematical homogenization theory, different kinds of coupled models have been investi- 
gated/derived. Besides the combination of fluid flow and convective-difFusive transport, the cou- 
pling among different kinds of species by chemical reactions have been discussed for example 



12| . see also the references cited therein. Further cross couplings of the water flow by heat, 



chemical or electrostatical transport are studied formally in 4] . It is worth pointing out a totally 
different context, where a nonlinear coupling quit e analogous to the one of our problem occurs - 
the phase-field models of AUen-Cahn type, see jlOj for more details on the modeling, analysis, and 
averaging of such models. Investigations concerning variable scaling and their influence on the limit 

. . . . ri 

equations is illustrated (by means of formal two-scale asymptotic homogenization) in [3], where 
different choices of ranges of the Peclet number are considered. In the same spirit, but this time 
rigorously, different scale ranges are examined for a linear diffusion-reaction system with interfa- 
cial exchange in (24|. Moreover, hybrid mixture theory has been applied to swelling porous media 



with charged particles in ^ and 



Formal upscaling attempts of the Nernst-Planck-Poisson sys 



2l|. It is 



tem using formal asymptotic expansion are reported, for instance, in [J|, |17| . [20| and 
worth pointing out that [20] and 2l| succeed to compute (again formally) microstructure effects 
on the deforming, swelling clay. In spite of such a good formal asymptotic understanding of the 
situation, rigorous homogenization results seem to be lacking. Only recently, Schmuck published a 
paper concerning the rigorous upscaling of a non-scaled Stokes- Nernst-Planck-Poisson system with 
transmission conditions for the electrostatic potential, ^28^. Furthermore, Allaire et al. studied the 
stationary and linearized case in Q. Our paper contributes in this direction since we perform the 
rigorous homogenization of the SNPP system for different boundary conditions as well as for vari- 
able choices of scalings in e, where e is a scale parameter referring to a (periodically-distributed) 
microstructure. The main focus of the paper thereby lies on the investigation of the influence of the 
boundary condition and scalings in e on the structure of the effective limit equations. This paper 
is built on 25]. However, we corrected essential errors concerning the use of Poicare's inequality. 
Furthermore, we introduce suitable redefinitions of the electrostatic potential in order to provide 
a more clearly arranged form of our homogenization results. Most important for the applications, 
we extend our results for different choices of boundary conditions for the electrostatic potential 
and include Stokes equations to our analysis in order to describe the interactions with the fiuid 
fiow. 

The paper is organized in the following way: In Section [31 we present the underlying micro- 
scopic model equations - the Stokes-Nernst-Planck-Poisson system. This is the starting point of 
our investigations. The Nernst-Planck equations describe the transport (diffusion, convection and 
electrostatic drift) of and reaction between (number) densities of colloidal particles. The electro- 
static potential is given as a solution of Poisson's equation with the charge density which is created 
by the colloidal particles as forcing term. The fluid flow is determined by a modified Stokes equa- 
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tion. Basic results concerning existence and uniqueness of weak solutions of this coupled system 
of partial differential equations are stated in Theorem 13.71 in Section [31 Moreover, Section [3] con- 
tains the definition of the basic heterogenous and periodic geometric setting. The (small) scale 
parameter e introduced here balances different physical terms in the system of partial differential 
equations and plays a crucial role in the homogenization procedure. Furthermore, e independent 
a priori estimates are shown for both Neumann and Dirichlet boundary conditions of the electro- 
static potential in Theorem 13.51 and Theorem 13.61 In Section |4l we state the basic definitions and 
well known compactness results concerning the method of two-scale convergence. The main idea 
is to obtain an "equivalent" system of partial differential equations that can reasonably describe 
the effective macroscopic behavior of the considered phenomena. We achieve this by investigating 
rigorously the limit e — !■ using two-scale convergence. Our analysis focuses on the influence of 
the choice of the boundary condition for the electrostatic potential and the different choices of 
scalings in e on both the a priori estimates and the structure of the limit problems. The main 
calculations are included in Section [4. II and Section [4.2l The crucial point is the nonlinear coupling 
of the system of partial differential equations by means of the electrostatic potential, and therefore, 
the passage to the limit e — > in the nonlinear transport terms of the Nernst-Planck equations 
and the Stokes equation. The main result (Theorems 14.51 14. 7[ and Theorems I4.11[ 14.131 14.15P 
of the paper discuss for which choices of scaling we can pass rigorously to the limit e — > 0. The 
results of this homogenization procedure and the structure of the limit equations are emphasized 
in Remarks H El [7] and HI El E] and in Section O 

2. The Underlying Physical Model 

We list in Table [T] all variables and physical parameters that are used in the following including 
their dimensions. Thereby, L is a unit of length, T a unit of time, M stands for a unit of mass, C 
for a unit of charge, while K represents the unit of temperature. 

In this section, we formulate a system of partial differential equations describing colloid dy- 
namics. Following e.g. [ll| and [scl, we impose to our system the balance of mass as well as 
the conservation of electrostatical charges. Note that in most applications, colloidal particles are 
charged [s^. Besides standard transport mechanisms (convection and diffusion), a charged dis- 
persion of colloidal particles is also transported by the electrostatic field created by the particles 
themselves as well as by the possibly charged soil matrix. Further interaction potentials (e. g. van- 
der-Waals forces or an externally applied electrostatic field) may also act on the colloidal particles. 
Throughout this paper we neglect the latter effects and focus on the investigations of the intrin- 
sic electrostatic interaction. Following Chapter 3.3 in js^, the positively (-I-) and negatively (-) 
charged particles are modeled in an Eulerian approach by some number density , which is trans- 
ported by the total velocity that consists of two parts: First, the convective velocity term v^'^'^'^ 
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V 


[L/T] 


velocity 


p 


[M/L/T^] 


pressure 


V 


[M/L/T] 


kinematic viscosity of the fluid 


P 


[M/L^] 


density of the fluid 


c 


[1/L3] 

[LVT] 


number density 


D 


diffusivity 


V 


H 


outer unit normal 


<P 


[V] [ML^/T^/C] 


electrostatic potential 


a 


[ML/T'^/C] 


surface charge density 


z 


[-] 


charge number 


e 


[C] 


elementary charge 


eoCr 


[C/V/L] 
[ML^/T^/K] 


dielectrostatic permittivity • relative permittivity 


k 


Boltzmann constant 


T 


[K] 


absolute temperature 



Table 1: List of the variables and physical parameters and their dimensions. 



due to the fluid flow within the porous medium in which the colloidal particles are transported. 
This is the same for all types of charge carriers. Second, the drift term -ydrift,* ^ is different for 
both kinds of charge carriers, can be calculated from the drift force = _2:^eV<? via 



„drift,± 



± rdritt.i 



with proportionality coefficient / and an electrostatic interaction potential In applications, / 
is sometimes also called electrophoretic mobility and is related further to the diffusivity by 



the Stokes-Einstein relation 



fcT ' 



„± ^ „drift,± 



[30[. The total velocity can therefore be expressed by 



- V 



hydr 



kT 



Inserting this expression into the standard convection-diffusion-reaction equation for a number 
density results in a modified transport equation which is also known as Nernst-Planck equation. 
On the boundary F of the considered domain Q we assume no-flux condition, which supplements 
the so called "no penetration" model, described in Together with an appropriate choice of 
the initial conditions c^'*^, the transport of the charged particles can be described properly by the 
following equations: 



dtc^ 



kT 



kT 



c^V^j = i?* (c) 



in(0,T)xi7, (la) 
on(0,r)xr, (lb) 
in {t = 0} X n. (Ic) 



with c := (c+,c^). The right-hand side in the Nernst-Planck equation include chemical reac- 
tions between the particles, source terms et cetera. 

The electrostatic interaction potential has to be calculated using Poisson's equation ((2a|) . The 
effect on the electrostatic field implied by the charged particles themselves is included as right-hand 
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side. This equation may be supplemented by Neumann or Dirichlet boundary conditions which 
correspond to the surface charge and the so cahed C potential of the solid matrix, respectively. 
Depending on the application in the geosciences either of the boundary conditions is given for 
example by measurements. 

— (z+c+ - Z-C-) in (0, T) x (2a) 

W<P-i^^a on (0, T) x Tjv, (2b) 

<P = <Pd on(0,T)xrD. (2c) 

In order to determine the fluid velocity v^^'^'^ we solve the modified Stokes' equations for incom- 
pressible fluid flow ([5a|) . I3bl As force term on the right hand side we take into account the drift 
force density. These equations are supplemented by a no slip boundary condition. 

-TjAv'^y^' + -Vp - --(z+c+ - z-c-Wif''^ in (0, T) x Q (3a) 
P P 

V-v^y^'^O m{0,T)xQ (3b) 

^'^''^'■ = on(0,r)xr. (3c) 

Remark 1. (Part of) the system ([1]), ([2]), dS]) arises in more general contexts. It plays a role when 
determining ion distributions (for example around colloidal particles or in a ion channel) and also 
in the framework of semiconductor devices especially if the convective term is neglected. We refer 
the reader to [l^, [2^ for aspects on the modeling and analysis of the semiconductor equations. 

3. Pore Scale Model 

In this section, we incorporate the physical processes described in Section [2] in a multi-scale 
framework and state basic properties of weak solutions as well as results concerning solvability 
of our problem. On the one hand, the phenomena considered in Section [2] take place on the mi- 
croscale and, on the other hand, the physical behavior we are interested in occurs on a macroscopic 
domain. In the framework of colloids, the transport takes place within the pore space of a porous 
medium that is defined by its soil matrix. The definition of the idealized underlying geometry 
which characterizes the highly heterogenous porous structure is depicted in Figure [TJ The (small) 
scale parameter e is introduced to scale/balance the different terms in the governing system of 
partial differential equations ([TJ, ([2]) and ([3]). 

Let us consider a bounded and connected domain i? C M", n E N with an associated periodic 
microstructure defined by the unit cell Y — (0, 1)". In the following we only consider the physically 
meaningful space dimensions n £ {1, 2, 3}. The unit cell Y is made up of two open sets, see Figure[TJ 
The hquid part Yi and the sohd part Y^ such that Y; U = F and YiHYs = 0,1^ n ^ T. 
Especially, the solid part does not touch the boundary of the unit cell Y and therefore the fiuid 
part is connected. We call e < 1 the scale parameter and assume the macroscopic domain to be 
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e 



Figure 1: Standard unit cell (left) and periodic representation of a porous medium (right). 

covered by a regular mesh of size e consisting of e scaled and shifted cells F/ that are divided into 
an analogously scaled fluid part, solid part and boundary. Let us denote these by Yf^^ K/j, and 
r^^i, respectively. The fluid part/pore space, the solid part and the inner boundary of the porous 
medium are defined by 

f2, :=[]Yi% f2\n,:^[]Yl,, and :=|Jn,,. 

i i i 

Consequently, since we assume that f2 is completely covered by e-scaled unit cells Y^ and, in 
particular, since the solid part is not allowed to intersect the outer boundary, i. e. dQ ~ ^. 

The objective of the paper is to rigorously investigate the limit e 0. The focus thereby lies 
on the coupling between the colloidal transport, the fiuid fiow and the electrostatic potential. We 
weight the different terms in ([T|), ([5]) and ^ with the scale parameter e in order to derive reasonable 
macroscopic model equations. In the framework of colloids, a non-dimensionalization procedure 
which can be used to motivate the choice of scaling has been done for example in jsoj]. However, 
since the system ([T]), ([2]) and ^ is used to describe various kinds of applications, different choices of 
scaling may be interesting depending on the underlying physical problem. We focus on the infiuence 
of the nonlinear coupling of the SNPP system due to the electrostatic potential and therefore regard 
Neumann as well as Dirichlet boundary condition for the Poisson equation and consider only the 
scaling of the coupling terms. For the ease of presentation, we assume that D := Z?+ = 
and z := z+ = — and suppress here the (constant) parameters r], p, z, e, k,T, D, er, eo as well 
as the superscript ''^^^"^ within all the equations. The resulting system of scaled partial differential 
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equations is referred here as Problem P^: 



— e ZiWe + Vpe — 




m (^U, J j X J4, 


(4aj 


V ■ Vg- — 


u 


in i^u, -i j X J /£ , 




We = 


n 
u 


on (^u, J j X (^i £ u Oii) , 


(4cj 


-e°'A<^^ = 




in (0,T) X n,, 


(4d) 


e^y^e ■ V ^ 




on (0,T) X P. AT, 


(4e) 


= 




on (0,T) X r,^D, 


(4f) 


e°V<?e • V = 





on (0,T) X dn, 


(4g) 


Vc± T £^C±V<Pe) = 


Rt{ct,c-) 


in (0,T) X 12., 


(411) 


^ ± e'^c^V^'e) • V = 





on (0, T) X (/; U df2) , 


(4i) 


ct- 




in {< = 0} X n^. 


(4j) 


vity constraint — 


c~ = 1 which is 


quite standard for the system 





see e.g. [26l]. This constraint can be relaxed in the case of Neumann boundary condition for the 
electrostatic potential, for the homogenous case see also 27 [. 



Remark 2. We could add a variable scaling also for the convective, diffusive and reactive terms. 
However, we concentrate on the role of the electrostatic potential The same choice of scaling 
in the equations for cf is especially justified in the case that both types of particles have similar 
properties except of the sign of the charge. On the outer boundary df2 we assume homogenous flux 
conditions for the concentration fields and the electrostatic potential as well as no slip boundary 
conditions for the velocity field. However, different linear boundary conditions could be chosen 
instead without notable changes in the calculations. For a discussion on different boundary con- 
ditions on the inner boundary and their influence on the results of the homogenization procedure 
we refer to the discussions in Remark SI [5l [7] and [3 [U [TT] and in Section [Sj 

Multiplying the system of equations ^ with the test functions ipi £ (iJp (i?^))" , (/?2, </53, V' G 
H^{n^) and integrating by parts we get the following weak formulation of Problem P^: 

/ e^Vue • Vipi - peV ■ (fiidx = / -e^(c+ - c~)V<Pe ■ Vidx (5a) 

/ Ve ■ VV'dx = (5b) 

/ e°'V<l'e-Vip2dx- I £"V^e • l^(/J2 rfOa; = / (c+ - c^) (^2 dx, (5c) 

{dtcf,^s)iH^y.H^+ f {-v,cf +S/ct±e^cfV^,) -Vipsdx^ [ Rfict .cj)^zdx. (5d) 



Definition 3.1. We call {ve,Pe,^e,cf ) a weak solution of Problem ifv^ e L°° {{),T; Hq{Q^)^ , 
Pe e L°° (0, T; L^in^)) , <?>e e (0, T; H^f2^)) and cf e (O, T; L^{f2,)) n (O, T; H\f2,)) 
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with dtcf £ (0,T] {H^{Q^ jy^ and equations ^ are satisfied for all test functions (pi € 

We modify the drift term in the Nernst-Planck equation by replacing the coneentration fields cf 
with the cut off functions cf := max(0, cf): 

dtcf + V • {v,cf - Vcf T s^cfV^e) = Rfict,c;) in (0, T) x i7„ (6a) 

{-Vecf + Vcf ± e''cfV^^) • J/ = in (0, T) x {r, U df2) , (6b) 

cf = c±'° in {t = 0} X ^2^. (6c) 

The modified system consisting of ^ and (|^a|) -([^gl) is referred here as Problem P^. The weak 
solution of Problem is defined analogously to Definition 13.11 

Remark 3. The weak solution of Problem is also a weak solution of Problem P^. Furthermore, 
all non- negative weak solutions of Problem P^ are also weak solutions of Problem P^ . As stated 
in Theorem 13.71 Problem P^ has a unique solution which is the non-negative one. Therefore both 
problems are equivalent. 

To be able to state a result on the existence and uniqueness of weak solutions of Problem P^ , 
we assume the following additional restrictions for the ease of presentation. Especially item 2 
and 4 can be relaxed. Note that, e. g., nonlinear monotonic reaction terms can be handled using 
homogenization theory as treated in [l4j |. 

Assumption 1. 

1. On the geometry: We assume a perforated domain as introduced in Section\^ i. e. the pore 
space J7e is bounded, connected and has C'^'^ -boundary. 

2. On the rate coefficients: The reaction rates are assumed to have the following structure 
R^(cf ,c^) — ^{cf — c^). Especially, they are linear and employ conservation of mass for 
the concentration fields. 

3. On the initial data: We assume the initial data to be non-negative and bounded independently 
of e, i. e. 

< c^'°(x) < A for all x^fl. 
Furthermore we assume the following compatibility condition for the initial data, i.e. 

[ c+-^-c-'°dx= [ ado^ 

If a — this implies global electro neutrality for the initial concentrations. 

4. On the boundary data: We assume the boundary data a and <Po to be constant. 

In order to ensure unique weak solutions, we additionally require 
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Assumption 2. If the electrostatic potential <!> ^ is determined via the equations {jd!^ , |-^ep and J^gp , 

we assume the potential to have zero mean value, i. e. J^^ <Pg dx — 0. Furthermore, we assume 
the pressure to have zero mean value, i. e. J^^ p^ dx = 0. 

Theorem 3.2. Let {vi^,p^,(l>i,,cf ,c^) be a weak solution of Problem in the sense of Defini- 
tion \3.1\ Let furthermore Assumption[l\hold. Then the total mass M — cf + c'^ dx is conserved. 

Proof. We test the Nernst-Planck equations (jSdp with (^3 = 1, sum over ± and insert the 
structure of the reaction rates according to Assumption [1] which directly gives the statement of 
Theorem [221 

Theorem 3.3. Let {v^,p^,(!>^,c'^ ,c'^ ) be a weak solution of Problem P^. Let furthermore Assump- 
tionUihold. Then the concentration fields are non-negative, i.e. are bounded from below uniformly 
in e. 

Proof. We test the Nernst-Planck equations (jSdp with ip^ = (cf)- := min(0, c^) which yields 

/ dtcf{cf)_-v,cf-\7{cf)_+\7cf-\7{cf)_±e''cf\7<P,-\7{cf)_dx^ [ Rf{cf)_dx. 

The drift term cancels directly due to the definition of the cut off function . The velocity term 
cancels by standard calculations due to the incompressibility and no slip boundary condition. After 
summation over ±, we have 

= / ~{ct - c-){ct)- + ic+ - c-)(c-)- dx. 

We consider the reaction term Iji :— — (c+ — c^)(c+)_ + (c+ — c^)(c^)_ for the following cases: 

1. c+ > 0,c- > 0: /i? = 

2. c+ < 0, c- > 0: Ir = -(c+ - c-)c+ < 

3. c+ > 0,c7 < 0: Ir = (c+ - c-)c- < 

4. c+ < 0, c- <0:Ir^ -(c+ - c-)c+ + (c+ - c-)c- - -(c+ - c^)^ < 
In any case we have the estimate Ir < and therefore 

Ijt + llfe")-lli^(r..)) + (l|Vc+||i.(,,^) + ||Vc-||i.(,,^)) < 

Gronwall's lemma implies the statement of Theorem 13.31 since the initial concentrations are non- 
negative according to Assumption [TJ 

Theorem 3.4. Let (u^jp^, c^, c~) be a weak solution of Problem P^. Let furthermore Assump- 
tion]^ hold. Then the concentration fields are bounded from above uniformly in e. 
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Proof. The statement of Theorem 13 .41 follows directly from Theorem 13.31 combined with the vol- 
ume additivity constraint cf + = 1. The boundedness of the concentration fields can be 
proven in the case of Neumann boundary conditions for the electrostatic potential without the 



volume additivity constraint, for the homogenous case see Lemma 3.3.6. in [27| where Moser's iter- 
ation technique is applied formally. This formal proof can directly be extended to non-homogenous 
boundary conditions and linear reaction rates as defined in Assumption [T] and also be made rig- 
orous. A rigorous approach using Moser's iteration can be found for general nonlinear equation in 
[l^. However, an alternative and more straight forward way is to show that a maximum principle 
applies in the case of homogenous Neumann boundary conditions for the electrostatic potential. 
Since the solutions cf of Problem are non negative, cf can be replaced by cf in the the 
Nernst-Planck equations. Using ip^ = {cf — := max(0,c^ — A) as test function, we obtain 

/ dtcficf - A)+ ~ v,cf ■ V(c± - A)+ + Vc± • V(c± - A)+ ± e^cfWP, ■ V(c± - A)+ dx 

= f Rficf~A)+dx 

The velocity term cancels by standard calculations due to the incompressibility and no slip bound- 
ary condition and it remains 

±e^ I {cf - A)V^, ■ V{4 - ^)+ + AS/^e ■ V(c± - A)+ dx = [ Rficf - A)+ dx 

We consider the drift term separately. Using the identity [cf — yl)V^e • V(c^ — A)j^ — V<?e • 
^V(c^ — A)\ and integration by parts, leads to 

±6-^ / i(-Z\<P,)(c± - A)\ + ^e^(-Z\<?,)(c± - A)+ dx 

Here the homogenous Neumann boundary condition for the electrostatic potential prevents the 
occurrence of boundary terms. Summation over ± and inserting the Poisson equation leads to 

/ e^((c+ - c-)\{ct A)\ (c+ - c7)i(c- - A)\ + ^(c+ - c7)(c+ - A)^ 

—A{c^ — cj){cJ—A)^\dx=: / T^dx 

' J 

We now distinguish the following cases: 

1. c+ < A,c- < A.Td=0 

2. c+ >A,c-< A: Td = {c+ - c-)\{c+ ~ Af + A{c+ - c-){c+ - A) > 

3. c+ < A,c~ >ATd = -(c+ - c~)\{c- - Af - A{c^ - c-)(c- - /I) > 

4. c+ >A,c-> A 

Td = \{ct Af - l(c+ - A){c- Af \(,c- ^)(c+ - Af + i(c+ - .l)^ 

+ A[ct-c-f>^ 



10 



Here we used the identity (c+ — c^) — (c+ ~ A) — (c^ — A) and applied Young's inequality 
(3, 3/2) which leads to a cancelation of all but the last term. 

We now consider the reaction term := — (c+ — c^){cf — + (c+ — c^)(c^ — ^)+ for the 
following cases: 

1. c+ < A,c- <A:Tr = 

2. c+ > A, c- <A:Tr = -{c+ - c-)(c+ - yl) < 

3. c+ < A, c- >A:Tr = (c+ - c-){c- -A)<Q 

4. c+ > A,c- >A:Tr = -(c+ - c-)(c+ - A) + -(c+ - c7)(c7 - yl) = -(c+ - c,-)^ < 
Finally, since To > and Tr < we have 

i| (||(c+ - yl) + ||i.(„^) + ||(c7 - Vl)+||i.(,,^)) 

+ (||V(c+ - ^)+||i.(,,^) + ||V(c- - /l)+||i.(,,^)) < 

Gronwall's lemma implies the statement of Theorem 13.41 since the initial concentrations are 
bounded from above by A according to Assumption [TJ 

In the following Theorem we state a priori estimates that are valid if we assume Neumann 
boundary data for the electrostatic potential on Fi^. This corresponds to a physical problem in 
which the surface charge of the porous medium is prescribed. 

Theorem 3.5. Let Assumption]^ and\^ be valid. The following a priori estimates hold in the case 
of pure Neumann boundary conditions for the electrostatic potential: 

£"\\^e\\Lmo,T)>cn^} + e"||V<l'e||L2((o,T)xa) < C. (7) 
In the case /3 — a > 0, it holds 

ll«e||L2((o, T)xr2J + £||VVe||L2((o,T)xf2,) < C. (8) 

// additionally 7 — a > is fulfilled, it holds 

n^l^x J|c-||i2(r2^) + max ||c+||i2(f2,) + ||Vc7||i2((o,T)xr2.) + W"^ 4 hmo. T)xf2,) 

+ \\dtc^\\L^o,T;iH^{n,)y) + \\dtc^\\L^o.T;{HHn,)y) < C, 
In (0) and (0), C € ]R_|_ is a constant independent of e. 

Proof. To derive the a priori estimates we test ((5c)) with the potential (p^ which leads to 

e"llV<Pe|li2(r2,) < £|kl|L2(n)||<?e||L2(r.) + ||c+ - \\mn,)\\^e\\L^n^:) 

< V^MLHr,)C iW-PeWmn,) + s\\W<Pe\\mn,)) + ||c+ - cJ\\L2i^n^)\\W<Pe\\LHf2^) 

< C^(lklU^(r,) + 114 - c-|U2(,,^)) ||V^>,|U2(,,^). 
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Here we used e||<?e|||2(/^^) < C {\\^e\W2(Q^) + £^ll^^£|lL2(r2-)) '^^^^ some constant C independent 
of e, see [3| Lemma 3, Poincare's inequality for functions with zero mean value (cf. [5]) and £ < 1. 
This results in 

e'^W^'PeWmn.^ < C + ||c+ - c7|U.(r,,)) < C, 

since a is constant and the concentration fields cf are bounded uniformly in e, see Theorem 13.41 
Using once again Poincare's inequality leads directly to statement ([7]) after integration with respect 
to time. The constant C remains bounded e-independently due to Theorem l3.4l and Assumption [TJ 
We test (jSaP with the velocity field and apply Poincare's inequality for functions with zero 
boundary values, i.e. ||<(5e||2 < Cpe||V(^e||2 with some constant Cp independent of e, see 
page 52. This leads due to the incompressibility of and the e-independent boundedness of cf 
according to Theorem 13.41 to 



This results in 

e\\Vv,\\L2^n^^ < ef^C\\V<P,\\L2^n^) < C, 

if /3 — a > 0, since the right hand side is bounded independently of e due to the estimates derived 
for the electrostatic potential. Using once again Poincare's inequality leads directly to statement 
(|8]) after integration with respect to time and the constant C remains bounded e-independently. 

In Theorem 13.41 we have already shown that c+ and are bounded by A uniformly in e. We 
test the Nernst-Planck equation (jSdp with ip^ = cf to obtain an energy estimate. This allows to 
bound also the gradient of the concentration fields. 

<[ \e^cf\/^,-Vcf\dx+ [ Rfcfdx<Ae''\\\7^,\\L2(n.j\\Vcf\\L2in.j+[ Rfcfdx 
< e'^'^'^Cs (llalli.(p^) + ||c+ - c7||i.(,,^)) + <5|| Vc±||i.(,,^) + Rfcfdx 

Here we used the estimate for the electrostatic potential derived above and that the velocity term 
cancels due to incompressibility of the fluid and the no slip boundary condition and Young's 
inequality. Summation over ±, sorption with 5 < 1/2 and estimation of the reaction terms via 

-(c+ - c-)c+ + (c+ - c-)c- < -{cf - c-f <0 finally leads to 

^ ^ -|-m2 I ||„-||2 \ I ^ /||Y7„+I|2 I ||Y7„-||2 



L2(r2,) + Ike IIl2(j7,) j + 2 (ll^'^s 11^2(^2,) + 11^^^ llL2(f2,)^ 



Integration with respect to time gives an uniform estimate of the gradient if 7 — a > since a is 
constant and the concentration fields are bounded independently of e. 
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To conclude the proof of Theorem l3.5[ we still need to derive estimates for the time derivatives 
dtcf of the concentration fields. By the definition of the {H^)' norm and by equations (jSdp . we 
obtain 

\\dtct\\{H^{f2,)y ^ sup {dtcf,if}(^H^y,H' 

¥>eHi(f2J,||¥'llHi(r2,)<l 

< sup ((|lc+ - c;\\L2i^n,) + A\\ve - e^W<P,\\L2^n,-, + ||Vc±||i2(^2,)) MH^n,)) 

¥>eHi(^^=).ll'p!l„i(j2,)<i 

< Wcth^n^) + \\c;\\mn,) + /Lllt^ellL^cfi,) +^£''""e"l|V<Z',|U2(,,^) + ||Vc±|U2(„,) < C, 
if 7 — a > due to the uniform estimates for the gradient of the concentration and the potential 
derived above, respectively. Integration with respect to time therefore yields the last statement of 
Theorem 13.51 

In the following Theorem we state a priori estimates that are valid if we assume Dirichlet 
boundary data for the electrostatic potential on F^. This corresponds to a physical problem in 
which the surface potential of the porous medium is prescribed. In application in the geosciences 
this boundary condition is related to the specification of the so called C potential. We define 
the transformed electrostatic potential := <P^ — <Pd- Since is a constant according to 

Assumption [U ^hom f^ig^g ^j^q following set of equations: 

in(0,r)x/2e, (10a) 

in(o,r)xr„ (10b) 
in(o,r)xar2. (lOc) 

Theorem 3.6. Let Assumption{l\ be valid. The following a priori estimates hold in the case of 
Dirichlet boundary conditions on for the electrostatic potential 

e"-2||^-^^-||^2((o,T)xI7.) +e"-'l|V<^r'"||L^((0,T)xr2.) < C. (11) 

In the case ^9 — a + 1 > 0, it holds 

\\Ve\\Lmo,T)xO,) + e||VWe||L2((o,T)xf2,) < C. (12) 

In the case 7 — a + 1 > 0, it holds 

max Ijc7||i2(f2,) + max ||c+||i2(fi^) + || Vc7|li2((o,T)xr2.) + II Vc+|li2((o,T)xr2,) 

+ \\dtct\\L^a,T-{H^n,)y) + ll^tc^ 11^2(0^7^.(^1(^2^)),) < C. 
In and \13\) . C G M+ is a constant independent of e. 

Proof. We test equation (jlOap with the translated potential cl>^°™ and use Poincare's inequality 
for zero boundary data, see [12|. This leads to 

^"l|Vf^>°'"|li2(,,^) < l|c,+ - crlU2(,,^)||<?>^-°|lL^(r2.) < l|c,+ - c-|U2(,,^)£Cp||V^>^"|U2(^^), 



^hom ^ Q 
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which results in 

e"-i||V<?^--|U.(a) < Cp\\ct - c^Wmn,). < C 

Here we have used the boundedness of the concentration fields provided by Theorem 13.41 with 
C being a constant independent of e. Using again Poincare's inequality leads to 

— 2 1 1 ^hom 1 1 ^ 

Altogether, we obtain the statement ([TT]) directly after integration with respect to time. By means 
of Theorem l3.41 the constant C remains bounded e- independently. 

The rest of the statement in Theorem 13.61 follows analogously to the proof of Theorem 13.51 
since due to the definition of the translated electrostatic potential and Theorem 13.61 it holds 

e"-i||V<l>,|U2(„^) -e"-i||V<?hom||^,^^^^^ < c 



The (stationary) system consisting of (jlap and (j2ap without convective term is well known 
as drift-diffusion model or van-Roosbroeck system in the theory of semiconductor devices 
Analytical investigations treating existence and uniqueness of solutions of this system can be 



found in [19| and 



26[. Extensions of the svstem (|lal) and (Pal) to the Navier-Stokes equations have 



27[. The results proven there can be carried 



been considered analytically, for instance, in [26|, 
over to system ([4]) and the following Theorem holds true: 

Theorem 3.7. Let Assumptions^ and [7] be valid. For each e > there exists a unique weak 
solution of Problem in the sense of Definition \3.l[ 

4. Upscaling of Problem 

This section is the bulk of the paper. Here we pass rigorously to the limit e — > in the non- 
stationary pore scale model P^ for both the Neumann and Dirichlet case and different choices 
of scaling (a,/3,7). For this aim we apply the method of two-scale convergence which has been 
introduced by Nguetseng in [23] and further developed by Allaire in [l!]. An introduction to this 
topic and the application of this method to basic model equations can be found, for example, in 
Q and For the reader's convenience, we state the definition of two-scale convergence as well 
as^he basic compactness result for functions defined on a time-space cylinder, see, e.g., [3l and 

Definition 4.1. A sequence of functions {y^e} ((Oj'^) ^ ^) said to two-scale converqe to 

a limit (po belonging to ((0, T) x i? x Y) if, for any function il> in D ((0, T) x i7; (^)) , we 
have 

lim / / ip^{t,x)^] [t,x, — ] dxdt= I / ipo(t, x,y)ip(t, x,y) dy dx dt. 
''^oJoJo ^ £' JoJny.Y 
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2 

In short notation we write ipe —^ipo. 

A sequence of functions {^Pe\ in L'^ {{0,T) x F^) is said to two-scale converge to a limit (pQ belonging 
to i2 ((0, T)x f2x r) if for any function i/; in D ((0, T) x C^^^{r)) , we have 

i-T I- , . i-T 



lim £ 



93e(i, (i, a;, — ) dox dt = / / ipo{t, x,y)'tp{t, x,y) dy dx dt. 



Here D ((0, T) x I?; C^,.(r)) and D ((0, T) x C;j^,(r)) denote the function space of infinitely 
smooth functions having compact support in (0,r) x fl with values in the space of infinitely 
differentiable functions that are periodic in Y and respectively. The following compactness 
result allows to extract converging subsequences from bounded sequences and therefore yields the 
possibility to pass to the two-scale limit provided that suitable a priori estimates can be shown. 

Theorem 4.2. 

1. Let {(pe} be a bounded sequence in i^((0,T) x f2). Then there exists a function (po in 

((0,r) X fi X Y) such that, up to a subsequence, (p^ two-scale converges to ipo- 

2. Let {ipi;} be a bounded sequence in L^ (0,T; {f2)^ . Then there exist functions (po in 
L^ (Q,T;H^ and (pi in L^ ((0,T) x SI; Hp^^{Y)'j such that, up to a subsequence, (p^ 
two-scale converges to po and Vipe two-scale converges to Vx'i^o + Vj,(y9i. 

3. Let {(pe} and {e'Vip^} be bounded sequence in L^ ((0, T) x fi). Then there exists a function LpQ 
in L^ ((0, T) X i7; H^^^ (F)) such that, up to a subsequence, {pe and e'S/kpe two-scale converge 
to ipo and Vyipa, respectively. 

4. Let {ipe} be a bounded sequence in L^((0,r) x Tg). Then there exists a function ipQ in 
L^ ((0,T) X [2 X r) such that, up to a subsequence, (pe two-scale converges to ipo. 

Proof. For a proof of the time independent case we refer e. g. to [3| and 
easily be carried over to the time dependent case. 



23| . The proof can 



One difficulty is that the a priori estimates that have been derived in Theorem 13.51 and The- 
orem 13.61 are at first only valid within the perforated domain . Therefore an extension of the 
functions u^, Vue,Pe, V<^e, c^, 9tc^, Vc^ is necessary, such that appropriate a priori estimates 
can be extended and that the limits for e — >■ can be identified in function spaces on fl. This 
procedure is quite standard and we refer to Q, B, B> H and Q 

for the strategy and the proof 

of the following 

Theorem 4.3. For the concentration fields cf we apply a linear extensions operator E (z C {H^{Q^), 
such that 
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is valid. 

The pressure field p;, is extended via 



E{Pe) := 



Pe in i^e, 

/y-e Pedy in each Y^^, 
and the following uniform a priori estimate holds if we assume zero mean value in fl: 



l|£^(P£)l|L2((0,T)xf2) < C. 

The other variables are extended by zero into fi. Then fl^ can be replaced by f2 in the a priori 
estimates from Theorem \3.5\ and Theorem [ 



However, for the ease of presentation we suppress the notation of the extensions and write 
again (^^ instead of E{ip^). 

In the next two subsections we consider the homogenization of system Q for both the Neumann 
and Dirichlet case via two-scale convergence. The statements on the two-scale limits of the extended 
functions and on the derivation of the macroscopic limit equations are deduced using the a priori 
estimates in Theorem 13.51 and Theorem 13.61 Special attention is paid to the coupling via the 
electrostatic interaction and the influence of the ranges of scaling on the limit equations. We first 
state the following 

Definition 4.4. We define the averaged macroscopic permittivity and diffusion tensor by 

A, / {S^J + dy^^j (y)) dy, (14) 

where tpj are solutions of the following family of cell problems (j = 1^ . . . ,n) 

-Ay(pj (y) = inYi, (15a) 

Vyipj {y) ■ V = —Cj ■ V on F, (15b) 

ipj periodic in y. (15c) 

We define the averaged macroscopic permeability tensor by 

K,,^ I w]dy, (16) 

where Wj are solutions of the following family of cell problems (j — 1^ . . . ,n) 

—AyWj + VyiTj = Cj in Yi (l^a) 

Vy-Wj^O in n xYi (17b) 

Wj =0 inYs (17c) 

Wj periodic in y (17d) 
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Furthermore, we define the following cell problem 

-Ay^{y)^l inYi, (18a) 

^{y) = {) our, (18b) 

Lp periodic in y. (18c) 

^.1. Neumann boundary condition 

We define := £"<l>e- 
4-.1.1. Homogenized Limit Problems for Poisson's Equation 

Theorem 4.5. Let the a priori estimates of Theorem \ 3.5\ be valid. Then the following two-scale 
limits can be identified for the electrostatic potential and its gradient V<?e ; There exist functions 
3oeL'^{0,T;H^{Q)) and 3ieL'^{{0,T) x Q;Hp^^{Y)) such that, up to a subsequence, 

3,{t,x) ^ i'oit,x), 

VS,{t,x) ^VMt,x)+yy^iit,x,y). 

Proof. We consider the estimate ([7]) in Theorem 13.51 which imphes 

W^eWmn) + ||V<^e||L2(r2) < C. 

Theorem 14.21 ensures the existence of the two-scale limit functions. 



Theorem 4.6. Let {vg,pg,(l>i,,cf ,c'^) be a weak solution of Problem Pg in the sense of Defni- 
tion \3.1\ Assume that cf converge strongly to in L^ ((0,T) x f2). Then the two-scale limits of 
due to Theorem \4-5\ satisfy the following equations: 

-V, • (i?V,<^o(i,x)) - (To = \Yi\ {c+{t,x)^CQ{t,x)) m (0,T) x 

D\/^3o{t, x)-iy = on (0, T) x dH. 

Proof. To prove Theorem l4.6l we test Poisson's equation ([5c|) with test function (tpo{t, x)-^e^i (t, x, 
which leads to 

VS^(t,x) ■ V (^ijjo{t,x) -\-eipi (t,x, dxdt 



Jf2 

T 



£^ [i^oii-, + £ipi (t, X, — )) dxdt 



/ Xe {x) {cf {t, x) — cj {t, x)) I -00 {t,x) + eipi it,x, — j] dxdt 
We then pass to the two-scale limit e — !■ using the properties we have stated in Theorem 

i-T 



X{y) {Vx^o{t,x) +Vy^i(t,x,y)^ ■ {Vx'ipoit.x) +Vyipi{t,x,y)) dydxdt 



JnxYi 

T 

aoipo{t,x) doy dx dt 



J sixr 

/ X (y) {ca {t,x)- Cg {t,x)) ipo {t, x) dy dx dt 
J nxYi 
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Now, we choose 'ipo{t,x) — 0, which leads, after integration by parts with respect to y, to 

-Vy ■ (wMt,x) + S7y<l>i{t,x,y)'^ =0 in {0,T)x fix Yi, 

(wMt,x)+S/ySi{t,x,y)'^ -ly^O on {0,T) X X r, 

(l>i{t,x,y) periodic in y 

and, therefore, also to 

-AySi{t,x,y) ^0 m{0,T)xQxYi, (19a) 

VySi{t,x,y) ■ V = -V^io{t,x) ■ V on {0,T) X f2 X r, (19b) 

'l>i{t,x,y) periodic in y. (19c) 

Due to the linearity of the equation, we can deduce the following representation of 

^1 (i , X, y ) = ^ (y) d,^ 3o{t,x) (20) 
j 

with (fj being solutions of the standard family of j = 1, . . . , n cell problems ([T5|l. 

On the other hand, if we choose ipi{t, x, y) = 0, we may read off, after integration by parts with 

respect to a;, the strong formulation for • 

Vx-I^j Vx^f){t,x) + Vy^i{t,x,y)dy^ - j a^doy ^\Yi\{c'^{t,x) - CQ{t,x)) in(0,T)xl2, 

(^J ^^Mt,x)+WySi{t,x,y)dy^ -,^ = on (0,T) x dH. 

Inserting the representation ([20)1 of <l>i yields 

V:, • (^DV^i>o{t,x)^ -ao = \Yi\ (c+(t, x) - c^(t, x)) in 

DV ,.3o{t, x) ■ v = ondn 

with diffusion tensor D being defined in (jl4p and (Tq := Jp ctq doy. 

Remark 4 (Modeling of <?o). In the case a = 0, it follows = (l>^. Therefore, we have an 
macroscopic equation for the leading order potential which is directly coupled to the macro- 
scopic concentrations c^. The case a < implies that and V<?e converge to zero. However, for 
any a an effective equation can be derived for the limit <2>o of . 

4.I.2. Homogenized Limit Problems for Stokes' Equation 

Theorem 4.7. Let the a priori estimates of Theorem \S.5\ he valid, i.e. especially /3 > a. Then the 
following two-scale limits can be identified for the velocity field and the gradient eVwg.' There 
exists vo e ((0,T) x [2;H^^^{Y)) such that, up to a subsequence, 

Vs{t,x) vo{t,x,y), 
e\7ve{t, x) VyVo{t, X, y). 
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Proof. We consider the estimate (|8]) in Theorem 13.51 which impHes due to Theorem 14.21 the 
existence of the two-scale hmit functions. 

The convergence for are standard, see and we follow directly the procedure there in- 
cluding the right hand side which is due to the electrostatic interaction. Depending on the choice 
of the scale range, this possibly leads to a coupling of the flow with the electrostatic potential and 
the concentration fields as stated in the following 

Theorem 4.8. Let {v^,pi;,(l>g,cf ,cj) be a weak solution of Problem in the sense of Defini- 
tion \3. il Assume that cf- converge strongly to Cg in ((0,r) x Q). 

For (3 > a the two-scale limit of due to Theorem \4-T\ satisfies the following equations: 



. {4{t,x)-Co{t,x))V,Mt,^), /3 = «l\ . .^TWO 

j[t,x) ^ -K[\/xPo(t,x) + < } m[0,T)xf2, 

V [ 0, /3 > a J / 

V^-va{t,x) ^0 in{0,T)xf2. 

Proof. Choose eip if, x, |) as test function: 

eVv^{t, x) ■ e^VtJj (t,x, — j — x)eV ■ ip (t,x, — j dxdt 

£^^^{c^{t,x) - cj(t,x))V$^(t,x)ip (^t,x, ^) dxdt. 



Passage to the limit leads to 



Jn 



Poit, X, y)Vy ■ %l){t, X, y) dy dx dt = 0, 

J ny.Y 
which gives po{t, x, y) = po(i, x). 

We define the space = {Vy • = 0, ■ /y^ V c^J/ = Oj = on (0, T) x i? x Yg} and choose 
ij{t,x, |) e as test function: 

eVve(i, a;) • eS/ip (^t, x, — ^ — Peit, x)V • (t, x, — ^ dx 



Ja \ E' 



a 



[ -e'^{cj{t,x) - cj{t,x))y$e{t,x)i) (t,x,-) dx 

Passage to the limit leads to 

7 Vy.,.W,^-p.^..^dydxJ^^^^^^-~^''"^-''^^^^^^^ 
JQxYi I 0, P > a 

Here we applied that ip €V^, i.e. Wy-tp — Q holds. The property po = Pq{x) yields 

/ ~PQ{t,x)Vx ■ [ I 'tp{t,x,y)dy\ dxdt^O 
Jn \Jyi J 
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Integration by parts inserting the properties of the orthogonal of T/^ and identification of the 



pressure po as m 



12| leads to 



-AyVa + V^pa + VyPi = < > m (0, T) x 12 x F; 

[ 0, /3 > a J 

Vy(;o=0 in (0,r) X 12 X F/ 

Va:- vody = in (0,r) x 

vody-u^O on (0, T) x 9^2 

Wo = on (0, T)x Six K,. 

If /? = a, we define the modified pressure pi — pi + {cq —Cq)<1'i in order to determine a macroscopic 
extended Darcy's Law. Due to the linearity of the equations vo can be represented as 



Y, 



(t, x,y) ^ Wj{y) ( dj:^po{t, x) 



3 



{c+{t,x)-Cg{t,x))d.j:^SQ{t,x), 13^ a 
0, I3> a 



with Wj being solutions of the cell problems (fT7|). We define the averaged velocity field via 

vo{t,x)^ / vo{t,x,y)dy. (21) 
Jy, 

which leads, after integration with respect to to 

/ \ (ct(t,x) - c;;(t,x)W.j,So(t,x), I3^a\\ 

V [ 0, /3 > a J / 

V^-vn{t,x)^0 m{0,T)xQ 
with the permeability tensor K being defined in (jl6p . 

Remark 5 (Modeling of vq). In the case /3 = a, we derive an extended incompressible Darcy's 
law. Besides the pressure gradient, an additional forcing term occurs due to the electrostatic 
potential. In the case /? > a, the electrostatic potential has no influence on the macroscopic 
velocity, which is then determined by a standard Darcy's law. 

4.I.3. Homogenized Limit Problems for the Nernst-Planck Equations 

Theorem 4.9. Let the estimates of Theorem \3.5\ be valid. Then the following two-scale limits can 
be identified for the concentration fields cf and their gradients Vc^ in the case 7 — a > 0.' There 
exist functions Cg (t,a;) G {{0,T);H^ {f2)) and ci{t,x,y) G ((0,T) x [];H^^^{Y)) such that 
(up to a subsequence) 

cf{t,x) Cg 
Vcfit.x) ^V^c^{t,x)+Vycf{t,x,y). 
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Proof. The statement of strong convergence holds true due to the extension of the concentration 
fields with the properties defined in Theorein l4.3l and Aubin-Lions compact embedding lemma. 

Remark 6. The strong convergence of the concentrations cf in (0,T;L^ (J7)) enables us to 
pass to the limit e — ^ also in the convective and drift term of the Nernst-Planck equations (|5dp. 



Theorem 4.10. Let {v^,pi;,(l>g,ct ,cj) be a weak solution of Problem in the sense of Defini- 
tion \3.1\ Assume that V<Pe o.nd two-scale converge as stated in Theorem \4-5\ and Theorerr \4-7\ 
respectively. 

Then the two-scale limits of the concentrations as stated in Theorem satisfy the following 
macroscopic limit equations: 

Dc^{t,x)\'.^$o{t,x), -f^a 
0, 7 > a 

= \Yi\R^{c+it,x),c^it,x)) m {0,T)xf2, 

.Dct(t,x)Vx^o(t,x), 7 = a \ 
vo{t,x)c^{t,x)- D\/^c^{t,x)± { ov^ J ^ \ \ -1^ = on {0,T)xdQ, 

0, -y > a \ J 

Proof. We choose (p2,3 = ipo{t,x) + si^iitjX, as test function in the Nernst-Planck equa- 
tions ([5d| and obtain: 

I -c^(t,x)dt {il^o(t,x) -\~ eil;-t {t,x, ^ 



{t,x)dt (^ipo{t,x) +eipi (t,x, + {-Ve{t,x)cf{t,x) 



J n 

-\- \'cf{t,x)±e^cf{t,x)V<Peit,x)) ■ V (xpa{t,x) + eil^i (t,x, dxdt 

J Rf{cf{t,x),cj{t,x))(^ipo(t,x)-\-eipi(^t,x,—^^ dxdt. 

Due to Theorem 14.91 and Assumption [U we pass to the two-scale limit e — > 0. 

/ -c^{t,x)dtiJoit,x) + i ~ vo{t,x,y)c^{t,x) {yc^{t,x) -\~\7ycf{t,x,y)) 
JnxYi \ 

{Co (Va^'^o + Vj^^i), 7 = a 1 \ 
> ■ {'^ xi^G{t,x)+^ yipi{t,x,y)) dydxdt 
0, 7 > a J / 

{c'^{t, a;), Cq {t, x))'i(jo{t, x) dy dx dt. 

In the case 7 = a we define c^ :— ± c^Si. We choose ipo = 0, which leads, after integration by 
parts with respect to y, to 

-Ayct{t,x,y) = m{Q,T)xnxYi, 

j_ j_ \ c^(t,x)'S/x^oit,x) ■ ly, 7 = a 

Vycf{t,x,y)-,y = -VxC^{t,x)T{ ' \ on(0,T)xl?xr, 

[ 0, 7 > a J 

Ci{t,x,y) periodic in y. 
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Due to the hnearity of the equation, we deduce the fohowing representations for cf: 

ct{t,x,y)=Y,^,{y)dx^ct{t,x)±{ « ^ ' ^ K (22) 
J [ 0, 7>aJ 

where ipj is the solution of the standard cell problem ([T5|) . 

On the other hand, if we choose il>i{t,x,y) = 0, we read off the strong formulation for c^, after 
integration by parts with respect to x, and after inserting the representation (j22p of cf: 



\Yi\dt4{t,x)+Wx ■ [vo{t,x)4{t,x)-DWx4{t,x)± 



0, 7 > a 

= \Yi\Rt{c+{t,x),c^it,x)) in (0,T)x/?, 

i 4- I -Dcn Vj^^o, 7 = a I \ 

vo{t,x)c^{t,x) - DV.xC^{t,x)± { " }] -1^ = 011 {0,T)xdf2, 

[ 0, 7>«J/ 
with D and vq being defined in (fH)) and (HI]), respectively. 

Remark 7 (Modeling of Cq). Mainly two different types of limit equations arise for the macro- 
scopic problem description. In the case j = a, the transport of the concentrations is given by 
Nernst-Planck equations. Thereby the limit of the electrostatic potential and vq are given in 
Theorem 14.51 and Theorem 14.71 The upscaling procedure then yields a fully coupled system of 
partial differential equation. In the case 7 > a, the electrostatic potential has no direct influence 
on the macroscopic concentrations. The equations for the concentrations simplify to a convection- 
diffusion-reaction equation. Depending on the choice of /3, the effective equations might be coupled 
only in one direction. 

The two families of cell problems ([T5|) and (|15p yield the same solutions and therefore the same 
macroscopic coefficients (up to the constant parameters that we have suppressed for the ease of 
presentation). 

4-2. Dirichlet boundary condition 

4.2.1. Homogenized Limit Problems for Poisson's Equation 

We define "P^ :— e"-2(^hom -^^j-^jj^jj fulfills the following set of equations: 

-e^Aie = ct -cj in (0,T) x Q^, (23) 

^e^O on(0,r)xr„ (24) 

e'^vS^ -ly^O on (0, T) X df2. (25) 



Theorem 4.11. Let the a priori estimates of Theorem \3.b\ be valid. Then the following two-scale 
limits can be identified for the electrostatic potential and the gradient eV<?e ; There exists <l>o€ 
L^{{0,T) X f2;H^^^{Y)) such that, up to a subsequence, 

S,{t,x) ^ 3oit,x,y), 
eVS,{t,x) ^ VySo{t,x,y). 
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Proof. We consider the estimate ([7]) in Theorem 13.51 which imphes 

ll^e||L2(r2) +e||V<^e||L2(f2) < C. 

Theorem 14.21 then ensures the existence of the two-scale hmit functions. 



Theorem 4.12. Let [v^^p^^'P^^c^ ^c^) be a weak solution of Problem in the sense of Defini- 
tion \3.1[ Assume that cf converge strongly to Cq in ((0,r) x f2). Then the two-scale limit of 
<Ps due to Theorem satisfies the following equations: 



$o{t,x)= yj^^ipj{y)dyj (cJ(t,x)-Co {t,x)). 

Proof. To prove Theorem 14. 121 we choose (t,x, |) as test function in (^5]) which leads to 

/ / eVi>eit,x) -Ve^ [t,x,-\ dxdt ^ I I {cf{t,x)-cj{t,x))^(t,x,-)dxdt. 
Jo Jn \ eJ Jq Jq \ eJ 

We then pass to the two-scale limit £ — ^ using the properties we have stated in Theorem 14. Ill 

/ / [\7y^o{t,x,y) ■Vyilj{t,x,y)] dydxdt = / / (c^ {t,x) - {t,x)) ^{t,x) dy dx dt 
Jo JuxYi ^ 'Jo J nxYi 

After integration by parts with respect to y, the strong formulation for niay be read off: 

-Ay^fi{t,x,y) = c+{t,x) - c-{t,x) in {=,T) X n xYi, 

S = in (0, T) X i7 X r, 

<pQ periodic in y. 



Inserting the cell problem (|T8|) . we get 

^0 = J^^ <^o dy = (^J^^ ip dy^ (c+ - Cg ). 
Remark 8 (Modeling of ^o)- In the case a = 2, it follows = ^(p^-cpjj and therefore 

^0 = <?S°'" + '^D = J^^ + ^Ddy= (^j^^ pdy^ (c+ - Cq ) + \Yi\^d- 

The macroscopic representation is directly coupled to the macroscopic concentrations Cq. The 
case a < 1 implies that and V<?e converge to and zero, respectively. However, for any a an 
effective equation can be derived for the limit <Po of 'P^. 

4-2.2. Homogenized Limit Problems for Stokes' Equation 

Theorem 4.13. Let the a priori estimates of Theorem \3.6\ be valid, i.e. especially /3 > a — 1. 
Then the following two-scale limits can be identified for the velocity field and the gradient e ; 
There exists vo G ((0,r) x Q;Hl^^{Y)) such that, up to a subsequence, 

2 



V, 



\t,x) voit,x,y), 



2 

£Vwe(t, a;) ^yVo{t, X, y). 



23 



The convergence for is standard, see [12| and we follow directly the procedure there including 
the right hand side which is due to the electrostatic interaction. 

Theorem 4.14. Let (u^jp^, c+, c^) he a weak solution of Problem in the sense of Defini- 
tion \8.1[ Assume that cf converge strongly to in ((0,r) x i7). 

For j3 > a — 1 the two-scale limit of due to Theorem \4-13\ satisfies the following equations: 



vo{t,x)= / vo{t,x,y)dy = -KVxPoit^x) m(0,r)xX2, 

V^-vo{t,x)=0 m{Q,T)xQ. 

Proof. Choosing eip{x, |) as test function, it follows analogously to the proof of Theorem 
that Pq — Po{x) holds. 

Defining the space — {Vj, • ijj — 0, • / -0 dy = 0, V' = on (0, T) x i? x Yg} and choosing 
'(/'(x, -) e as test function, leads in the limit e to 



Iq If2 X - (4 - Co ) Vj;^>o V' rfa; dt, /3 = a - 1 



^ f 

JnxYi I 0, /3 > a - 1 



We now follow the proof of Theorem 14.81 Finally, integration by parts results in 
-AyVQit,x,y) + VxPo{t,x) +Vypi{t,x,y) 
In the case /3 = a — 1, we define the modified pressure pi = pi + (cj — Cq )^o. This allows to 



-{c^{t,x) - Cq {t,x))Wy'pQ{t,x,y), /3 = a - 1 
0, /3>a-l 



determine a standard incompressible Darcy's Law and finishes the proof of Theorem 14. 141 

Remark 9 (Modeling of vq). The fluid flow is determined by a standard Darcy's law. The is 
no direct coupling to the electrostatic potential, since it is only present in the modified pressure 
term pi . 

4-2.3. Homogenized Limit Problems for the Nernst-Planck Equations 

Theorem 4.15. Let the estimates of Theorem \3.6\ be valid. Then the following two-scale limits can 
be identified for the concentration fields and their gradients Vc^.' There exist functions Co{t, x) G 
((0, T); iJi (i7)) and ci{t, x, y) G ((q, T) x Q; H^^^ {¥)) such that (up to a subsequence) 

cf{t,x) co{t,x), 
Vcf(t, x) ^ VxCoit, x) + Vycf(t, X, y). 

Proof. The statement of strong convergence holds true due to the extension of the concentration 
fields cf- with the properties defined in Theorem l4.3l and Aubin-Lions compact embedding lemma. 

Remark 10. The strong convergence of the concentrations in L^ (O, T; L^ enables us to 
pass to the limit £ also in the convective and drift term of the Nernst-Planck equations (jSdp . 
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Theorem 4.16. Let (w^jp^, c^, c~) he a weak solution of Problem in the sense of Defi- 
nition \3.1\ Assume that V^g and two-scale converges as stated in Theorem \4-ll\ and Theo- 
rem \4-l^^ Then the two-scale limits of the concentrations as stated in Theorem \4-15\ satisfy the 
following macroscopic limit equations: 

|rz|9tc±(i,x) + V, • {vo{t,x)c^(t,x)-DW,c^{t,x)) = \Yi\R^{c+ {t,x),c^ {t,x)) m (0,r) x [2, 
{vo{t,x)c'^{t,x) - D\/^c^{t,x)) • 1/ = on (0,T) x dH. 

Proof. We choose ips — tpo{t,x) + eipi (t, a;, f) as test functions in the Nernst-Planck equa- 
tions ((5d)) . 

cf{t,x)dt (v^o(i,x) +£i/'i (^.2;,^)) + {v,{t,x)cf{t,x)+Wcf{t,x)±e'^cf'^^,{t.x)) 
• V {^o[t^ x) + ei/'i {t^ X, — ^ ^ dx dt 

J Rf{cf{t,x),c^{t,x))(^'ipo{t,x)+eipi(^t,x,^^^ dxdt. 

Passage to the Hmit e — ^ yields 

n-co (i, x)dtil^a {t, x) + {-va{t, x, y)cj {t, x) + (VcJ {t, x)-\-S/ycf {t, x, y)) 
_ -2xY 



Jf2 



± 



CQVy(Po{t,x,y)) ■ {V xipo{t, x) yijjiit, X, y))dydxdt, 7 = a - 1 

0, 7 > a - 1 



(cq {t,x), Cq {t, x) )'0o {t, x) dy dx dt. 

/o JnxY 



We define 



-± _ ± , J Co Vy<?o, 7 = a - 1 



0, 7 > a - 1 

and choose V'o = 0, which leads, after integration by parts with respect to y to: 

-Aydf{t,x,y) = in {0,T) X nxYi, 

Vycf {t,x,y)-v = - V^^Cq {t, x) ■ V on (0, T) x Q x P, 

c^{t,x,y) periodic in 

The linearity of the equation yields ()22p as representations for c^ supplemented by the family of 
cell problems (jisp . 

On the other hand, if we choose ^i(t, x, y) — 0, we read off the strong formulation for after 
integration by parts with respect to x and after inserting the representation (|22p of c^: 

\Yi\dtC^{t,x) + V., ■ ivoit,x)c^it,x)-DVxC^it,x) = \Yi\R^{c+it,x),c„it,x)) in (0,r) x f2, 
{~vo{t,x)c^{t,x)+DW.xC^{t,x)) -ly^Q on (0,T) x dH, 

with vo,D being defined in ((T4|) and ((2T|) . respectively. 
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Remark 11 (Modeling of Cq). The transport of the concentrations is determined by a convection- 
diffusion-reaction equation. The is no direct couphng to the electrostatic potential, since it is only 
present in the modified higher order concentration term ci. Depending on the choice of /3, the 
effective equations might be coupled only in one direction. 

5. Discussion 

We wish to point out the following aspects: 
In SectionlU we considered the rigorous passage to the two-scale limit e — > for different boundary 
conditions of the electrostatic potential and different ranges of the scale parameter (a, /?, 7) and 
have derived the corresponding two-scale limits of Problem P^. We classified conceptually difFerent 
types of limit systems. In all cases, auxiliary cell problems need to be solved to be able to provide 
closed-form expressions for the effective macroscopic coefficients. Depending on chosen model, the 
macroscopic problem is coupled only in one direction or fully coupled. Solving these problems 
numerically is computationally challenging due to the mass balances that have to be fulfilled and 
the diverse boundary conditions, especially periodic ones. The most crucial point is that an ap- 
propriate fixed point iteration has to be constructed depending on the nature of the nonlinear 
couplings. Moreover, corrector estimates will be needed in order to make it possible to compare 
the effective solutions/problem descriptions with the oscillatory solution/microscopic model. The 
different structures of the resulting effective equations of the homogenization process are under- 
lined in Remark 21 Remark [5l Remark [7] for Neumann boundary conditions for the electrostatic 
potential (i. e. given surface charge) and in RemarkUl Remark|9l RemarkfTDfor Dirichlet boundary 
conditions for the electrostatic potential (i.e. given potential). In the colloid literature, one can 
also find the so called perfect sink boundary condition for the concentration fields instead of the 
no-penetration boundary condition, i.e. = on (0, T)x F^. In the framework of homogenization 
this would lead together with the strong convergence of the concentration fields to = as limit. 
Obviously, this does not provide a suitable model for colloidal transport phenomena. 

The following question arises naturally: Given a particular scenario of colloidal transport in 
the soil, which is the best/most reasonable mathematical (limit) model that should be considered? 
Answering this question is not limited to choosing the precise values for the choice of the appropri- 
ate boundary conditions and the scale range (a,/3,7). It also requires a careful calibration of the 
model by an intensive numerical testing of the chosen set of limit equations. Further adjustment 
by experimental measurements and parameter identification procedure may need to be done to 
make the model quantitatively. 

It is worth noting that, using two-scale convergence, we could not pass to the limit e — > for all 
choices of the parameter ranges. However, in these cases formal two-scale asymptotic expansions 
can be applied in order to pass formally to the limit e — > using the transformation := 
exp (=F^) which arises especially when treating drift diffusion problems (compare, e.g. [3;Q|)- 
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An alternative is to treat a linearized system as has been considered via rigorous homogenization 
in the stationary case in 
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